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ABSTRACT 


Several free radical addition polymerizations exhibit the diffusion limited gel or 
Trommsdorff effect. The models available till recently predict the behaviour of batch reac- 
tors reasonably well but are inapplicable to semibatch operations or to reactions under non- 
isothermal conditions, operations which are of industrial importance. The recent model of 
Ray et al. (23) focuses on such reactors, but needs to be confirmed experimentally. An 
experimental set-up to study polymerization in such reactors has been assembled. Methyl 
methacrylate (MMA) is polymerized in a computer controlled, 1-liter Parr reactor. A se- 
ries of experiments on bulk polymerization of MMA under different temperature histories 
[near-isothermal, step decrease and step increase] at a fixed initiator (AIBN) concentration 
of 25.8 mol/m^, have been conducted. Two important process output variables, namely, 
monomer conversion and average molecular weight, as determined from intrinsic viscosity, 
have been measured at different sampling times during the course of polymerization. The 
conversion and molecular weight data are compared with the model predictions obtained 
using the experimental (non-isothermal) temperature histories. Reasonable agreement is 
observed between theory and experiment. It may be added that the parameters of the 
model of Ray et al. were evaluated using isothermal batch reactor data from the literature, 
and no additional curve fitting has been performed in the present study. 
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CHAPTER 1 


INTRODUCTION 


Ever since Norrish and Smith (1) and TrommsdorfF et al. (2) reported the presence of 
the gel or auto acceleration effect in free radical polymerizations, scientists have been trying 
to obtain a fundamental understanding of this phenomenon. The gel effect is characterized 
by a sudden increase in the rate of reaction as the conversion increases beyond a certain 
point. This effect is associated with an increase in the viscosity of the reaction mixture, 
which causes a decrease in the termination rate constant kt (= ktc + ktj) (see Table 1) 
due to diffusional limitations. This leads to higher monomer conversions than predicted by 
conventional models. As the reaction proceeds, the glass effect is also exhibited at higher 
conversions, when the glass transition temperature, Tj, of the reaction mass becomes 
higher than the reaction temperature, T. This decreases the propagation rate constant, kp, 
causing the reaction to stop short of complete conversion. 

Mathematical modeling of the gel and glass effects has been attempted by several 
workers (3-5) in the past few years. They modeled these effects using time-varying rate 
constants for the termination and propagation reactions. The bulk polymerization data on 
methyl methacrylate (MMA) published by Balke and Hamielec (6) allowed the tuning of 
the parameters in these models. Friis and Hamielec (3) and Ross and Laurence (4) used 
simple empirical correlations to describe the polymerization kinetics in terms of monomer 
conversion, temperature and free volume, and curve-fitted their parameters. Cardenas 
and O’Driscoll (5), on the other hand, considered the polymer species as a set of two 
populations, one having a degree of polymerization (DP) below some critical value and the 
other, having higher values of DP. Two different values of the termination rate constant, 
ki, were taken. These were curve-fit parameters based on previous experimental data (6). 
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All these models introduce an artificial break-point in either the monomer conversion or 
the chain length, or both, at which the diffusional restrictions become dominant. These 
and several other of the earlier theories have been reviewed by O’Driscoll (7) and Hamielec 
( 8 ). 

Diffusional limitation was introduced as an integral part of the termination and propa- 
gation reactions by Chiu et al. (9). The effects of composition, temperature and molecular 
weight were accounted for continuously. The phenomenological nature of this model led 
to several interesting studies (10-18), with reference to engineering applications, focusing 
on aspects such as parametric sensitivity (15), optimization (13,16), optimal parameter 
estimation (17,18), etc. More recently, Achillas and Kiparissides (19,20), extended the 
model of Chiu et al. using the diffusion theory of Vrentas and Duda (21), and the theory 
of excess chain-end mobility (22), to give a molecular basis for the theory. 

The model developed by Soong and coworkers for hatch polymerizations, though ap- 
plied to the simulation and optimization of semibatch reactors (13,14), is questionable in 
its extension to such operations. This doubt arises from the use of the initial value of the 
initiator concentration, [I]o, in the model to account for the molecular weight dependence 
of the diffusivity of polymer radicals. Thus, it is unclear as to what should be used in 
place of [I]o in the model after the sudden addition of initiator, monomer or solvent to 
the reaction mass. This is an important question, in view of the fact that in industrial 
reactors such additions are made routinely to obtain optimal operation. A similar doubt is 
present in the theory of Adiilias and Kiparissides (19,20), in which an empirical curve-fit 
parameter, jco, related, in turn, to the initial number average chain length, is in- 
troduced to give proper curve-fit of isothermal experimental data. Some work has been 
reported (18,23) recently which attempts to resolve this controversy. A new model has 
been developed (23) for the gel and glass effects whidi is applicable to both batch and 
semibatch reactors under a variety of operating conditions. The validity of this theoretical 
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model, however, needs to be established using experimental data under nonisothermal and 
semibatch reactor conditions. The present study focuses on the nonisothermal operation 
of batch reactors. 

Work on isothermal bulk polymerization of MMA has been reported earlier (6). These 
results are based on experimental studies on ampoule reactors in which mixing effects as 
well as heat and mass transfer limitations do not play an important role. In the present 
work, however, a 1 liter stainless steel Parr reactor (Parr Instrument Company, Moline, 
IL, USA) has been interfaced with a PC-XT and methyl methacrylate is polymerized in 
it under controlled temperature histories. Several physico-chemical phenomena which are 
present in industrial reactors (e.g. vaporization, heat transfer through the viscous reaction 
mixture, etc.) are also present in our experimental reactor, and so our experimental results 
can be used more meaningfully to develop techniques, models and design-methodologies for 
larger prototypes. Experimental data on the monomer conversions and average molecular 
weights have been obtained for a variety of (nonisothermal) temperature histories, such as 
step increase and step decrease, and the results compared against model (23) predictions. 

The implementation of a desired temperature history in polymer reactors is not a 
trivial task. Achieving good temperature control is not easy due to the unfavorable com- 
bination of high heats of reaction and very high viscosities of the polymerizing mixture. 
Added to this is the complication introduced by the unsteady nature of the batch reactor. 
Very little work has been reported on this aspect. Elicabe and Meira (24) reviewed the 
literature in this area, but a detailed discussion is not being presented here for the sake of 
brevity. Amrehn (25) discussed a simple model-predictive controller qualitatively, which 
manipulates the jacket liquid temperature to control the temperature of a polymerization 
reactor. More recently, the performance of several kinds of controllers on polymerization 
reactors hats been simulated by Congalidis et al. (26), and by Hidalgo and Brosilow (27). 
Soroush and Kravaris (28) have recently tested a globally linearizing control (GLC) algo- 
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rithm experimentally on a 5 liter stirred glass batch reactor in which MMA is polymerized 
under an optimal temperature history. They used the model of Chiu et al. (9) to compute 
this optimal history. Their control algorithm was found to track this temperature history 
very well. In our study, however, we decided to employ a simple proportional-integral (PI) 
algorithm with set-point changes, to control the reactor temperature in the desired man- 
ner. Actually, facility for both heating and cooling the reaction mass was provided, the 
rates of which could be manipulated using two different PI controllers. An empirical gain 
scheduling was also included in the control algorithm. This scheme was found to give a 
fairly rapid heating-up of the reaction mass when required, as for example, during the ini- 
tial period, or in experiments involving a step increase in the temperature at some point of 
time. The controller could also implement reasonably fast decreases in temperature, when 
so required. In addition, this controller could maintain the temperature constant whenever 
this was desired. However, the simple PI controller does not work too well when the rate of 
release of the (exothermic) heat of reaction becomes very rapid, some time after the onset 
of the gel effect, and simultaneously mixing becomes difficult due to the high viscosity of 
the reaction mixture. In spite of this defeciency, we were able to obtain experimental data 
uptil sufficiently high conversions and were able to validate the theoretical model of Ray 
et al.(23). 

In this thesis. Chapter 2 describes the experimental set up used along with details of 
computer interfacing to record and control the reactor temperature. Chapter 3 elucidates 
the dynamics of the control variables, the temperature controller design and the software 
to implement the controller. Chapter 4 gives the procedure followed to carry out the 
polymerization and the analytical methods adopted to estimate conversion and molecular 
weights. Results and discussions are presented in Chapter 5. Conclusions drawn from the 
present study and recommendations for future work are included in Chapter 6. 



CHAPTER 2 


EXPERIMENTAL SET-UP 


In this chapter the experimental set-up, which includes the Parr reactor assembly, the 
computer and its peripherals, sensors and their calibration, transducers and actuators, is 
described. 

2.1 REACTOR 

The polymerization was carried out in a one liter stainless steel Parr reactor (Fig.la) 
(Parr Instrument Company, Moline, IL, USA). This reactor is designed to operate at 
pressures upto 1.4 x lO^kPa and tempertures upto 350‘’C. It has split ring closures for easy 
access to the interior of the bomb without using heavy screw caps or cumbersome cover 
clamps. In this design the closing force between the head and the bomb is developed by 
six cap screws which are held in a hardened steel clamping ring that has been split into 
two halves. These rings slide into place from the sides of the bomb so as not to disturb 
any fittings attached to the bomb head. A steel retaining ring raised into place from the 
bottom of the bomb completes the assembly. 

The reactor assembly can be lowered into a furnace (Fig. lb) which is provided 
with electric heating coils (1.5 kW), surrounded by a thick layer of insulating material. A 
helical cooling coil is placed inside the bomb (Fig. Ic). The contents of the reactor are 
mixed with two turbine impellers each with six blades, placed one above the other on a 
stirrer shaft. The reactor is fitted with three valves for different purposes (Fig. la). Valve 
Vi is a gas release valve which can also be used to apply vacuum to the reactor. Valve 
V 2 is the ga 5 charging valve. Valve V 3 , connected in the same line as Vj, is meant for 
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collecting licjuid samples from the reactor when the latter is pressurized. Valves V2 and 
V3 are installed on a dip tube which extends to the bottom of the reactor. There is a 
pressure gauge (0 — 1.4 x 10'' kPa) connected to the reactor through a needle valve. This 
valve helps in isolating the pressure gauge while the reactor is being evacuated. There is a 
thermowell fitted to the reactor head which extends into the liquid section of the reactor. 
A non-conducting and thermally stable liquid, silicone oil, is used as the thermowell liquid. 
An Argon gas (IOL-AR2, Indian Oxygen Limited, New Delhi, India) cylinder is used to 
pressurize the reactor and to carry out the reactions in an inert atmosphere. 

A 1.25cm SS needle valve operated with a stepper motor (Unique Systems and Con- 
trol Private Ltd., New Delhi, India) is fitted to a cooling water line (1.25cm steel pipe) 
connected, in turn, to an overhead water tank, to enable control of cooling water flow rate. 
The flow rate of the cooling water is measured with an orificemeter (orifice diameter 0.3176 
cm). 

2.2 COMPUTER SYSTEM 

The system consists of an IBM compatible (WIPRO) PC-XT with a 640 kB RAM, 20 
MB hajd disk, one floppy disk drive, a DT-2805 I/O interface card (Data Translation Inc., 
Marlborough, MA, U.S.A) with its companion DT-707T screw terminal, and a printer. 


2.2.1 I/O INTERFACE CARD 

A DT-2805 low level analog and digital I/O interface card has been used for interfacing 
the process with the computer (Fig.2). The card provides 8 differential ADC (Analog to 
Digital Converter) channels with programmable gains of 1, 10, 100 and 500 and two DAC 
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(Digital to Analog Converter) ports. The two ports (ports 0 and 1 ) of the digital I/O 
section have 8 bits each for data transfer. These ports can be addressed separately or in 
combination as a single port by choosing a proper binary value (mask value). The analog 
channels are used to measure analog inputs like thermocouple signal, pressure transducer 
signal, etc. The digital ports are used to run a stepper motor and to switch a solid state 
relay. 

The I/O card can be programmed to perform 12-bit [resolution 1/(2^^ ~ 1)] a-nalog to 
digital (A/D) conversion at 40 kHz and digital to analog (D/A) conversion at 33 kHz. As 
this is a low level card capable of accepting signals as low as 20 mV (full scale), differential 
inputs become necessary to eliminate noise. 

All the outputs, analog or digital, are held constant until changed or the board is 
reset by the in-built zero order hold circuit. It has its own on-board clock and processor 
for fast operation. An external clock can also be used. The jumper settings have to be 
chosen carefully for different configurations. 

The interface card is mounted in one of the slots of the mother board of the PC. It 
draws 1.2 A at 5 V from the mother board. The DT-707T screw terminal is connected 
to the board with a ribbon cable. The sensor and actuator connections are conveniently 
made to the screw terminal which resides outside the computer. This screw terminal also 
supports the measurements for a set of standard thermocouples. Room temperature is 
measured with a built-in solid state sensor connected to channel 0 of the screw terminal 
and this temperature is used for making cold junction compensation. 
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2.3 SENSORS AND TRANSDUCERS 

2.3.1 TEMPERATURE MEASUREMENT 

There are four temperature sensors used in this experiment (Fig.2), which measure 
the reactor inside temperature, the cooling water outlet temperature, the reactor outside 
wall temperature and the cooling water inlet temperature. The first three of them are 
copper-constantan thermocouples, and the fourth is a solid state sensor connected with a 
temperature transmitter (Analog Devices Inc., Norwood, MA, USA, model 2B57A) and 
energized with a 28 V power supply. This solid state sensor has a span of 100°C in the 
range of 0 -140° C and was used to measure the inlet temperature of the cooling water. 
The 4-20 mA output of the sensor is converted into a voltage signal by connecting a resistor 
(500 ohm) in series in the circuit. The circuit diagram is given in Fig.3. 

The procedure for calibrating the solid state sensor is as follows. First, the range of 
the sensor is adjusted to be 0-100°C by trimming the zero and span of the analog device 
over the full range of voltage mea.surement (for example, 0- lOV) by dipping it in melting 
ice and in boiling water. Then, the voltage drop is measured as a function of temperature 
by immersing the sensor in a constant temperature bath (Julabo Labortechnik, Seelbach, 
Switzerland, model FlOMH) maintained at different temperatures, and measuring the volt- 
age after it stabilizes. The calibration information is shown in Fig.4. The data points were 
correlated using the following linear relationship: 

T{°C) = 13.9313V -27.8533 (1) 

where V is in volts. 

The temperature measurement by thermocouples is done by the subroutine ’ Mea- 
sureThermocouple ’. The details of the measurement procedure can be found in the PCLAB 
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manual (Data Translation Inc., Marlborough, MA, U.S.A) (29). Although the interface 
card supports the copper-constantan thermocouples, it was found necessary to recalibrate 
the individual thermocouples from time to time. This was achieved using the Julabo 
FlOMH constant temperature bath. The temperature read by the computer was cali- 
brated against the actual (bath) temperature and these data were fitted with a straight 
line and used for more accurate measurement. The least squares fit equations are; 

For the reactor temperature, 

Tactua( = 0.996114 r,nea*ure<f+ 1.08907 ' (2) 

For the cooling water outlet temperature. 

Tactual = 1.05097 Tmeasured + 2.96102 (3) 

For the reactor outside wall temperature 

Tactual = 0.99976 Treasured + 1.19875 (4) 

■These calibrations were checked from time to time. The above equations represent the 
temperature corrections made for the near isothermal 50°C run (NI50). This pseudo cali- 
bration information for the two thermocouples is shown in Figs.5 and 6. The thermocouple 
which measures the reactor temperature was calibrated along with the thermowell. This 
procedure allowed taking into account the static lag of the thermowell. The noise problem 
experienced with the thermocouples was greatly reduced after grounding the reactor with 
the ’ground’ of the I/O card. 
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The solid state temperature sensor had an excellent stability and reproducibility but 
its response was sluggish when compared to that of the thermocouples. It was therefore 
used to measure the inlet temperature of the cooling water whose variation was expected 
to be the least. A thermocouple was used to measure the outlet temperature of the cooling 
water. Even though the inlet and outlet temperatures of the cooling water were not used 
for control at present, these are available for future use in advanced control applications 


2 .3.2 FLOW MEASUREMENT 


The cooling water flow rate is one of the two manipulated variables for the control of 
the reactor temperature. An orifice meter was designed and fabricated in-house (Fig 7 ) to 
measure the cooling water flow rate. The pressure drop across the meter was measured with 
a diflferential pressure transmitter (DPT, Taylor Instruments, Faridabad, India) having a 
range of 0 - 10 cm. of water. The 4-20 mA current output of the DPT was converted into 
a voltage drop by passing it through a 500 ohm resistor. The DPT was energized with a 
24 V DC power supply. The circuit diagram of the DP cell is shown in Fig. 8 . 


For calibration, the range of the sensor was first adjusted to be between the minimum 
flow rate (4 mA) and the maximum flow rate (20 mA) by trimming the zero and span 
adjustments repeatedly until further adjustments were not required. Then the voltage drop 
was measured and recorded as a function of water flow rate in the system by changing the 
valve opening. The flow rate was found by measuring the water collected for a known 
interval of time. The flow rate (fr) vs. voltage drop (V in volts) information is given in 
Fig-9. The experimental data were fitted with two best-fitting polynomials in two zones 
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Fig. 7 Details of the orifice plate flow meter used for 

cooling water flow rate measurement 
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They are given by: 

For 2 < V < 2.7 

Jr{glmin) = 227.062F^ - 1946.25^^ + 5672.84\/ - 5342.18 (a) 

For 2.7 < V < 10.0 

fr{glmin) = -5.4522y^ + 152.036^ - 110.311 (6) (5) 

Where V is in Volts. 

2.4 ACTUATORS 

2.4.1 STEPPER MOTOR 


The cooling water flow rate was controlled with a needle valve (1.25 cm) actuated by a 
stepper motor (Unique Systems and Control Private Ltd., New Delhi, India). The shaft of 
the motor is rigidly fastened to the stem of the needle valve through a coupling. Since one 
revolution of the motor shaft is divided into 200 steps, each step rotates the stem through 
1.8° angle. One sequence of pulses will rotate it in the clockwise direction (closing the 
valve) while the reverse sequence will rotate it in the anticlockwise direction (opening the 
valve). These sequences, along with the binary values if connected to the first four bits of 
a port are shown in Tables 2 and 3. 

The TTL (Transistor Transistor Logic) output from the computer cannot run the 
stepper motor directly. The pulses from the computer are sent to a board having diodes and 
solid state relays connected as shown in Fig.lO which converts the pulses into a synchronized 



Table - 2 


Sequence of pulses to stepper motor for clockwise motion of the shaft 


Step No. 
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II 

III 

IV 
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12 

2 

0 

1 

1 

0 

6 
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Table - 3 


Sequence of pulses to stepper motor for anticlockwise motion of the shaft 


Step No. 

I 

II 

III 

IV 

Binary 

1 
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0 
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12 
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0 
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Fig. 10 Driving circuit for the stepper motor 
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pulse train. The stepper motor requires a 12 V DC power supply. The stepper motor driver 
program is given in Appendix-2. 

2.5 FINAL CONTROL ELEMENTS 

2.5.1 NEEDLE VALVE 

A mult-turn 1.25 cm needle valve Weis used as the final control element to control 
the flow of the cooling water. The stem of the valve is connected to the stepper motor 
by a rigid coupling which eliminates any backlash. The motor was mounted on a flexible 
cushion seat to allow up-and-down motion, while at the same time it was restricted from 
rotation around its axis by fixing it rigidly at its corners with long bolts. It takes nearly 
8 turns for the needle valve to open completely. As mentioned earlier, 200 steps of the 
stepper motor lead to one rotation of the stem. Therefore, 1600 steps are required to open 
the valve completely. Such a large number of discrete positions of the valve make the 
control action almost continuous. 

The stem position (in steps of the stepper motor, S) of the needle valve was also 
calibrated against the cooling water flow rate [fr(g/min )]. The experimental data were • 
fitted with the following two polynomials: 

For 15 < fr < 500 

S = 8.49779 X IQ-^fr^ - 6.76177 x IQ-* fr^ + 0.42872/r + 0.843042 (a) 

For 500 < fr < 900 

S = 6.74481 X 10““/r® - 1.59919 x lO-^rV 1.45021 X 20~V^^ 


-0.0615486/r2 12.1968 fr - 800.288 (b) 


( 6 ) 
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These equations were used for controlling the flow rate. Fig. 11 gives the information on 
calibration. 


2.5.2 FURNACE CONTROL 


The heater power input is time-proportioned and implemented using a solid state 
relay (JDA 331000, Electronic Relays, Bangalore, India). The DC terminals of the relay 
are directly connected to a digital bit of the DT-707T. This relay can be set ON (high) 
or OFF (low) by the interface board. In the AC terminals of the relay, the furnace coil is 
connected in series with the 230 V AC power mains regulated with a variac (see Fig.l2). 
A neon lamp was connected in parallel with the furnace which glowed to indicate the ON 
position of the heater. 




Fig. 12 Furnace control circuit 








CHAPTER 3 


REACTOR TEMPERATURE CONTROL 


The problems often encountered in industrial liquid phase bulk polymerization reac- 
tors are associated with heat released by exothermic reactions and increase in viscosity. 
The heat of reaction ranges from 40 to 85 kJ/mol causing an adiabatic temperature rise 
of about 200-400°C. Large generation of heat, coupled with low thermal diffusivities of 
the system often cause thermal runaway as well as significant spatial gradients. Since 
higher temperatures give products of low DP (degree of polymerization), this can lead to 
disastrous consequences in terms of product quality. Because of these additional features, 
polymerization reactors pose a challenging control problem. 

Efficient control of temperature during exothermic heat release calls for a powerful 
cooling system. An equally effective strategy is necessary for rapid heating of the system 
from room temperature to the required temperatures (50-90°C), and to compensate for heat 
losses. Therefore, a single control variable, namely, reactor temperature, can be controlled 
best in such cases by suitably manipulating two variables, viz., cooling water flow rate and 
heater power input. In the present set-up it was decided to control one of the two variables 
at any time, depending on the sign of the error, T,p-T, instead of lumping the two variables 
into one (28). 

3.1 OPEN LOOP PROCESS DYNAMICS 


The process response curves (30) were obtained for the two manipulated variables 
by conducting step tests to study the heating and cooling dynamics of the present system 
individually. These open loop dynamic tests were conducted using water in the reactor. 
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The process output of the step response test, reactor temperature, was fitted in the 

least squares sense to a first order-with-dead time model represented by: 

^ /t^pexp(-tjs) 

Xis) TpS + 1 

where 

Kp = process gain = -r c ^an^ in output temperature 

y ^ o change tn heater power/ coolant flow input V ' 

y(5), X(s) = process output and input in transformed coordinates (c) 

td = process dead time or transportation lag (d) 

Tp = time constant of the process (e) (7) 

The model parameters were obtained by using the Gauss-Newton technique with the 
Levenberg-Marquardt modification [code in Cuthbert (31) as modified by Maiti (32)]. 



3.1.1 HEATING CYCLE DYNAMICS 


The heater power input weis controlled by the computer using time proportioning. 
The manipulated variable is the fraction of time for which the heater (150 V AC from the 
variac) is ON during a cycle time (20 s in this case). In the long run, this is effectively 
equal to the control on the average power to the heater. Thus, for example, to give 50% 
average power to the heater, full power is supplied for 50% of the cycle time (for 10 s in a 
20 s cycle). 

To obtain the heating dynamics, the reactor was first brought to a steady state 
temperature by continued heating for 2 s in a 20 s sampling time. During this test no 
cooling water was circulated through the cooling coil. Steady state was achieved in about 
4 hr. Then, a step of 2 s /20 s was applied by changing the heating time to 4 s/20 s. The 
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response of the system was followed by recording the reactor temperature at every sampling 
instant. This was continued for about 3 hr. The system parameters were obtained as (33), 

Kp = 361.082°C' /full power (i.e. 20 s heating / 20 s cycle) 

Tp = 213.87 min 

= 7.556 min (8) 

where the superscript, /i, indicates heating. 


3.1.2 COOLING CYCLE DYNAMICS 

The reactor was first brought to a steady state temperature of about 70°C. Both 
heating and cooling loops were then turned off. The cooling water flow was then suddenly 
turned on at a rate of 200 g/min. The falling temperature of the reactor was noted at 
every 10s. This gives the required dynamics of the cooling process alone. Here again the 
response was modeled using a first order-with-dead time transfer function. The calculated 
parameters of the cooling dynamics are (33), 

K; = -0.2009°C / {g/min) 

Tp = 6.8251 min 

= 0.6756 min (9) 

where the superscript, c, indicates cooling. 

3.2 CONTROL ALGORITHM 

The main thrust in the present work was to study the reaction kinetics of MMA 
polymerization under neax-isothermal conditions and with step changes in temperature. 
For this reason, relatively simple but efficient control algorithms were used in this work. 
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The control law studied was a simple coordination law between heating and cooling. 
This gave fairly satisfactory results. The decision to heat or cool is taken depending on 
whether the reactor temperature is below the set point or above it, i.e., if e(i) [=T„-T(i)], 
the error signal at time t,- is greater than a number e, then the reactor is heated using the 
heating PI controller. On the other hand, if e(i) is less than — e, then the cooling water 
loop is actuated using the cooling PI controller. Here, e is a small tolerance (=0.25°C) 
to account for the noise in the temperature signal. If thelerroij is within e, then neither 
heating nor cooling is done and the controller waits for the next sampling time. 

The simple control strategy discussed above gave satisfactory performance in terms 
of rapid set-point tracking. The performance curves for two experimental runs with water 
alone and with polyvinyl alcohol solution in the reactor are shown in Fig. 13, using the 
value of the controller parameters as given in Table 4. As seen in the figures, the control 
scheme performed satisfactorily. A combined control variable like the net heat added to 
the system could have been used (28) instead, but was not attempted. 

It was observed during trial runs with water in the reactor that once the temperature 
in the reactor dropped, it took a considerable amount of time for it to reach the set point 
again. This problem is to be expected in view of the delay time for the heater being an 
order of magnitude larger than that for the cooling system. To overcome this, power was 
supplied to the heater for 0.5 s / 20 s during both the cooling and the waiting periods. 
The exact amount of heating was tuned in such a way that it compensated for most of the 
heat loss to the surroundings. 

A sampling time of 20 s for heating and 10 s for cooling was used in this study 
since the cooling dynamics was relatively faster. The derivative action was not included 
for digital control as the temperature measurements were relatively noisy and derivative 
action without an effective filter would have been useless. 
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Fig. 13 Performance of the tuned \ i* 
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Table - 4 Cohen-Coon and Tuned controller parameters 


Parameter 

Cohen- Coon 
values 

Tuned 

values 

i<j 

0.0708 

0.07 


23.435 

900.0 


-45.7 

-70.0 

n 

1.864 

90.0 


— 

1.0 

^cas 

— 

1.8 
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Another faster PI controller operates in cascade in the flow-loop to minimize the effect 
of pressure fluctuations in the water line (Fig.14). This controller operates with a sample 
time of 5 s (smaller sampling times worsened the control). This controller manipulates the 
stem position of the valve till the desired and observed flow rates match. The controller 
parameters were tuned on-line, and the values of and t/“'’ are given in Table 4. 

S.3 DIGITAL CONTROLLER 

All the controller algorithms were implemented using the velocity form of the control law 
because the position form of the algorithm requires initial values of the controller output 
which are not known in practice. Also in this form, one does not compute the absolute 
value of the control variable at any given instant but only the increment to the control 
signal over the previous value (30). This method protects the controller from integral 
”wind-up” or saturation. The algorithm for a PI controller is discussed below. 


3.3.1 PROPORTIONAL-^INTEGRAL CONTROLLER 


In the continuous time domain the algorithm for a PI controller is given by, 

m(t) = Kc e(i) -f — f e{t)dt (10) 

Ti Jo 

where m is the manipulated variable (furnace power / cooling water flow). Discretising 
Eq.lO at the and (z — 1)*^ sampling instants, and taking the difference of the two, one 
obtains 

Am- = iir,(e.- - e,-a) -f — e.-r 

■■■ n 


(11) 
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Block diagram for flow 


cascade control 
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where T is the sampling time. Eq.ll represents the PI controller in its velocity form. 

3.4 COHEN-COON CONTROLLER DESIGN 


The first estimates of the controller parameters were obtained using the Cohen-Coon 
relations (30) . These parameters were subsequently tuned on-line. The controller param- 
eters were further fine-tuned for rapid heating (to attain the desired set point faster) and 
cooling dynamics, using viscous solutions of PVA (Jams Chemicals, Bombay, India) in wa- 
ter (concentrations chosen somewhat arbitrarily, to simulate approximately, the conditions 
of the reaction mass in the gel effect region) by giving step changes in temperature until 
smooth dynamics were obtained. For the cooling cycle, the calculated parameters required 
only very slight adjustments. The Cohen-Coon settings and on-line tuned parameters are 
given in Table 4. Further tuning of the controller parameters was done from time to time 
to account for changes in the ambient temperature and inlet temperature of cooling water 
with changes in season. 

3.5 SOFTWARE DEVELOPED 

In this study the real-time software was developed in Turbo Pascal. The interfacing 
was done by a TURBO PASCAL driver using DT 2805 interfacing card and DT 707T 
screw terminal panel [PCLAB manual. Data Translation Inc., Marlborough, MA, U.S.A, 
(29)]. The program was developed in a modular form with many procedures meant for 
specific purposes. The program can be used for both dynamics and control. The flow 
control procedure closely follows the work of Jana (34). The flow chart and listing of the 
program developed are included in Appendix 2. Some of the important procedures are 
discussed below. 
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READJN-CONTROL : This procedure reads in the set point, mode of control (P, PI 
or PID), and the controller parameters for the heating, cooling and the flow cascade loops. 

READ-IN.^TEP ; This procedure reads in the desired values of the control variables, 
namely, heating time per cycle time, and cooling water flow rate during open loop testing. 
Once the steady state is reached with the given conditions, a step change in one or both 
of the manipulated variables can be given. 

DISPLAY -CONTROL : This procedure displays the temperature, control action (heat- 
ing time and/or cooling water flow rate) at each instant. 

DISPLAY -STEP : This procedure gives the position of the manipulated variables and 
the temperature at each instant when the step test is being carried out. 

FLOW-CONTROL : Once the cooling water flow rate is specified by the primary 
temperature controller, this procedure takes up the job of obtaining the desired flow rate 
with the aid of the cascade PI controller. 

ORIFICEMETER : This routine reads in the voltage drop from the channel connected 
to the DPT and calculates the flow rate from the calibration information stored in the form 
of a polynomial. 

ORFW-CTRL-ACTION ; This procedure manipulates cascade controller’s output 
with PI action in order to obtain the desired flow rate. 

VALVE-MANIPULATION : This routine calculates the stem position required for a 
given flow rate and operates the stepper motor accordingly. 
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S_GFILTER: This procedure smoothens the flow measurement data using seven mea- 
surements taken at intervals of 2 s. The Savitzky-Golay formula (35) is used. 

ON : This sets the bit 1 of port 0 high (1) to switch on the furnace through the solid 
state relay. 

OFF : This sets the bit to low (0) to switch off the furnace. 

The general purpose routines which have been developed are described below. 

DECO. INC : This "include” file has routines to read the system time interacting with 
DOS and a procedure to draw boxes with the given corners. 

STEP.INC ; This program runs the stepper motor in a given direction for the specified 
number of pulses. This program issues the next set of pulses after analyzing the previously 
existing digital value. This program is currently designed for running the stepper motor in 
0, 1, 2 and 3 bits of one of the ports. By suitably changing the mask and the set of digital 
value, any set of bits can be used. 

The program can be interrupted at any time. The next action can be chosen from 
the given set of options. There is provision to change the set points and parameters in 
between. Information about interfacing routines like MeasureThermocouple, OutputDigi- 
tal Value, Measure Volts, etc., are available in the PCLAB manual [Data Translation Inc., 
Marlborough, MA, U.S.A, (29)]. 



CHAPTER 4 


POLYMERIZATION PROCEDURE AND 
CHARACTERIZATION TECHNIQUES 


In this chapter the techniques of monomer and initiator purification and the procedure 
for batch polymerization in a Parr reactor are described. The analytical procedures used 
for determining the monomer conversion and the molecular weight of the polymer produced 
are also discussed. 

4.1 PURIFICATION OF CHEMICALS 

4.1.1 MONOMER 


LR grade methyl methacrylate (MM A) (Burgoyne, Burbidges & Co., Bombay, India) 
was purified prior to use (6,36,37). The monomer was washed three times with equal 
volumes of a 5% NaOH (GPR, BDH, Bombay, India) solution in water to remove the 
phenolic stabilizers [hydroquinone(HQ) or methyl ether of HQ] present in it. The mixture 
was allowed to settle in a separating funnel forming two layers, the bottom inorganic layer 
was discarded leaving the monomer. The MMA was then washed at least three times with 
equal volumes of distilled water to remove traces of NaOH. Settling times of at least two 
hours after each NaOH and water wash were allowed. The residual water was removed 
from the monomer by passing it through beds (placed in series) of molecular sieves (type 
30-541, 4-5°A, Linde Division, Union Carbide, USA) and silica gel blue (LR, coarse, mesh 
3-8, NICE, Kochi, Kerala, India). Fresh batches of regenerated silica gel and molecular 
sieves (by washing and drying in an oven at 70°C for 5 hrs) were used for every 250 ml 
of monomer. The regenerated drying agents (silica gel regains its original blue color after 
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drying) were stored in air tight containers. 

The monomer obtained after the above steps was distilled under vacuum (P=50 mm 
Hg) at about 30°C (room temperature). Ice cold water was circulated in the condenser 
(Fig.15). The distilled monomer was stored in a refrigerator over fused CaC^ (Loba Chemie 
Pvt. Ltd., Bombay, India) in a desiccator. The final yield of the pure monomer was about 
60% of the original sample. 


4.1.2 INITIATOR 


The initiator, azobis isobutyronitrile (AIBN) (LR, Ranbaxy Labs. Ltd., S.A.S Nagar, 
Punjab, India) was recrystaJlized from methanol (LR, S.D.Fine Chemicals, India) (6). A 
saturated solution of AIBN was prepared (at room temperature) in a test tube by mixing 
an excess of AIBN with about 30 ml of methanol. The solution was filtered (No.T ordinary 
filter paper) into another test tube. The filtrate was then cooled in a refrigerator to 
crystallize the AIBN. The crystals were recovered by filtration. The crystals were dried in 
a vacuum oven at room temperature (heating causes thermal decomposition). 

4.2 BATCH POLYMERIZATION 

A 500 ml solution of required initiator concentration in monomer was prepared in 
a one liter round bottom flask. The solution was thoroughly mixed to ensure complete 
dissolution of AIBN. 

The solution was transferred to the Parr reactor. The reactor was carefully closed 
with the head and sealed with the half-rings and the steel retaining ring. All the valves were 




Fig. 15 Vaccum distillation set-iip 
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closed. The reactor contents were degassed by applying vacuum (P= 5 mm Hg) through 
the gas release valve for about 30 min. The vacuum pump (model 150D, Leo Engineering, 
Kalletumkara, Kerala, India) was connected to the reactor through ice-traps to remove any 
condensables in the gas. After degassing, the reaction mass was pressurized with argon 
(IOL-AR-2, Indian Oxygen Ltd., New Delhi, India) charged through the gas-charge valve. 
After about 20 seconds of the commencement of this operation the valve connected to the 
pressure gauge line was opened. Once the desired pressure of about 900 kPa was achieved, 
the gas-charge valve was closed. Argon was re-charged during a run if the pressure fell 
below 850 kPa. 

All the peripherals connections (temperature sensors, etc.) to the reactor-PC set-up 
(Fig.2) were checked and activated. The controller parameters and desired set point were 
supplied as inputs to the control program. Samples were withdrawn at required intervals of 
time through the valve V3 (see Fig. la). The sampling times were decided by the expected 
dynamics of the system. Samples of about 3 ml were collected in 50 ml sample bottles 
containing weighed amounts of solvent (10 ml of benzene, Extra pure, S.D Fine Chemicals 
Ltd., Boisar, India) and trace amounts of inhibitor (MEHQ). The sample bottles were 
kept ice-cooled in order to freeze the reaction as soon as the sample was collected. Before 
drawing the sample, about 2 ml of the sample coming out of the port was discarded, this 
being the estimated hold-up of material in the sample line. This procedure ensured that 
the samples drawn were really representative of the reaction mixture. The copper tube 
connected to the sample collecting port was cleaned with commercial grade acetone and 
dried after each sample before reuse. Sample bottles were reweighed to calculate the mass 
of the samples drawn. The occurrence of the gel effect can be inferred from the sudden 
change in the consistency of the material flowing out. During this period, samples were 
collected as often as possible. It was found that the argon pressure of about 10 atm above 
the liquid reaction mass was unable to push out the liquid through the sample port shortly 
after the onset of the gel effect due to the very high viscosity. 
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The reaction was stopped after the lapse of the required time or as soon as the 
stirrer stopped due to the high viscosity of the reaction mixture. After the reaction was 
complete, the reactor pressure was released slowly by opening the gas release valve and 
the reactor was then opened for cleaning up immediately. Most of the material inside was 
removed mechanically. The residual polymer was dissolved using dichloromethane (LR, 
NICE, Kochi,Kerala, India). In some cases (specially bulk polymerization at 70°C) it was 
quite difficult to open up the reactor after the reaction was complete. The contents had 
to be heated to higher temperatures (130-140°C) and kept at that point for a suflSciently 
long time (30 min ) to melt the polymer near the walls of the reactor. This enabled pulling 
out the reactor head with the cooling coils, stirrer and thermowell along with most of the 
polymer, in which these had become embedded. A considerable amount of polymer was 
first removed mechanically, and then the remaining polymer was dissolved with solvent 
(CR) acetone followed by (LR) dichloromethane. 

4.3 ANALYSIS OF MONOMER CONVERSION 

The monomer conversion was estimated gravimetrically (6). The polymer was pre- 
cipitated from the solution in the sample tube by adding about 100 ml of methanol (non- 
solvent). The precipitated polymer was filtered through a dried and preweighed filter paper 
(Whatman 41, Whatman International limited, Maidstone, UK). A few drops of the fil- 
trate were tested with about 10 ml of methanol to ensure that complete precipitation had 
taken place. The polymer was dried at 50°C under vacuum (P=5 mm Hg) for about 24 
hr (tested to be sufficiently long to ensure complete drying). The mass of the precipitated 
polymer was found by weighing the filter paper again. The polymer samples were stored 
in a desiccator till molecular weight analysis was done. 
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4.4 DILUTE SOLUTION VISCOMETRY 


The average molecular weights, M„, of polymer samples were determined by dilute 
solution viscometry (38-40). About 0.25 g of the dried polymer was dissolved in 30 ml of 
filtered extra pure benzene. The solution was kept in a stoppered standard flask for more 
than 24 hr (in a refrigerator) to ensure complete dissolution. The solution was carefully 
examined for complete dissolution before use. Keen observation was required as the poly- 
mer particles are transparent, making their presence difficult to detect. This solution was 
filtered through a sintered-glass filter (G4, For 4, Vensil Glass Works, Bangalore, India) to 
remove minute suspended particles. Benzene was filtered for further solution preparation, 
for dilution and rinsing of the viscometer in a separate filter. The glass filters were cleaned 
with chromic acid and distilled water, and dried in an oven. 

15 ml of the polymer solution was pipetted into the viscometer (Ubbelohde viscome- 
ter No: 501 01, Schott- Ger ate, Hofheim, West Germany) (Fig.16). The viscometer was 
carefully mounted on a stand built in-home and placed in a constant temperature bath 
(maintained at 30°C, using an Indotherm and a PTIOO sensor, IT401 D, PI controller, 
Indotherm Instruments Ltd., New Bombay, India). Efflux time measurements were carried 
out as described below. Suction is applied to tube 1 (Fig.16) by closing tube 2 with a 
rubber stopper to fill the working capillary (7), time bulb and upper reservoir (9). Then 
suction is discontinued and tube 2 is opened again. This causes the liquid column to sep- 
arate at the lower end of the working capillary (7). The efflux times are measured for 
descent from marks Mj to M 2 . These flow time measurements were repeated at least four 
times, with a reproduceability within 0.5 s. The original solution was diluted by pipetting 
5 ml of benzene into the viscometer. To help complete mixing of the solution with the 
added solvent, air was bubbled through the mixture in the viscometer using the pipette 
filler. The next efflux time measurement was taken after about 40 min. This was to ensure 
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attainment of homogeneity of the sample as well as temperature equilibration. 

As further dilution could not be made inside the viscometer (capacity of only 20 ml), 
the solution was taken out. 10 ml of this solution and 5 ml of filtered benzene were pipetted 
back into the viscometer. This was mixed as before and placed in the constant temperature 
bath and efflux time measurements were carried out. Flow times were measured for one 
more dilution by adding 5 ml of benzene. In this manner four efflux times were taken for 
any polymer sample. After thoroughly cleaning the viscometer, efflux time for the pure 
solvent was also measured. 

Appendix 3 gives the details of the calculations performed to obtain the intrinsic- 
viscosity average molecular weight from the efflux time measurements and concentration. 
The listing of the Fortran program used for double extrapolation and one specimen plot is 
also given in Appendix 3. 



CHAPTER 5 


RESULTS AND DISCUSSION 


The main thrust of this work wets to provide an experimental verification of the kinetic 
model developed by Ray et al. (23) for polymerization of MMA under step changes in 
temperature. Such operations are often encountered in industrial reactors. 

Before conducting the actual experimental runs, we attempted to find out if tempera- 
ture gradients existed within the Parr reactor during operation. A relatively concentrated 
solution of polyvinyl alcohol in water, having a consistency which was about the same (by 
visual inspection) as the polymerization reaction mass in the gel-effect region, was taken in 
the open reactor and heated with stirring. The thermocouple in the thermowell was moved 
around and the temperature was recorded at several different positions. These are given 
in Table 5. It is observed that the temperature is relatively uniform within the reactor, 
with a mean of 50.56° C and a standard deviation of 0.39, for a set point of 50°C. 

A series of experiments on bulk polymerization of MMA with different temperature 
histories [near-isothermal (NI), step decrease (SD) and step increase (SI)] at a fixed initia- 
tor (AIBN) loading, of 25.8 mol/m^, have been conducted. Two important process output 
variables, namely, monomer conversion, x, and average molecular weight, M„, as deter- 
mined from intrinsic viscosity, have been measured at different sampling times during the 
course of polymerization. Table 6 lists the controller parameters in all the runs conducted, 
while the detailed results of the monomer conversion and molecular weights are given in 
Appendix 1. It may be added that the controller gain was changed on-line to obtain faster 
heating (or cooling) during step changes in the temperature. 


The first two experimental runs were conducted under near isothermal conditions at 



Table - 5 


Radial temperature distribution using PVA solution 


Set point Tip= 50°C 


Distance (radial) 

Location of 

Height from top of reactor (cm) 

from wall (cm) 

stirrer (cm) 

8 (cm) 

4 (cm) 



Temperature (°C) 

Temperature (°C) 

0 (wall) 

at center 

51.01 

51.04 

2 

at center 

50.81 

50.73 

4 

at center 

50.63 

50.48 

6 

at one side 

50.76 

49.96 

8 

at center 

49.99 

50.09 

10 

at center 

50.08 

50.36 

other (wall) 

at center 

50.98 

50.99 




Table - 7 

Controller parameters used for each run 


Sno. 

Run No. 

Kc' 




Comments 

1 

NI50a 

0.07 

900 

-70 

90 


2 

NISOb 

0.07 

900 

-70 

90 

- 

3 

NI70 

0.08 

900 

-70 

90 

KJ 



0.07 




changed at t=8 min 

4 

SDla 

0.08 

900 

-70 

90 


5 

SDlb 

0.07 

900 

-70 

90 

Kl 



0.08 




changed at t=30 min 

6 

SD2 

0.07 

900 

-70 

90 

KJ, K' 



0.08 


-60 


changed at t=30 min 



0.07 




changed at t=80 min 

7 

SI 

0.07 

900 

-70 

90 




0.08 




changed at t=120 min 


CO.* 


Kc =1 
rr=1.8 
for all runs 
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50°C and TO^C (Run Nos. NISO and NI70, respectively). The conditions are refered to as 
near isothermal due to the fact that the reactor system takes about 6 to 8 min to attain 
the desired polymerization temperature (Figs. 17a and 18a). Vaporization of the monomer 
is prevented by a applying pressure of about 1400 kPa absolute, using an inert gas such as 
argon. The argon pressure also helps in driving out samples of the reaction mass through 
a port. 


The monomer conversion (z) for the two NISO runs (NISO a and b) are shown in 
Fig. 17b. The average molecular weight (M„) based on intrinsic viscosity measurements are 
shown as a function of conversion for the NISO b run in Fig. 17c. The duplicate runs were 
made to establish reproducibility of the results. Since the rise time (time taken to reach 
50°C from room temperature) for the NISO runs is very small, the conversion vs. time results 
are comparable to those obtained by Balke and Hamielec (6) on ampoule reactors (in which 
mixing and heat and mass transfer effects are negligible) under isothermal conditions. The 
excellent agreement between our conversion data (x vs. t) with those of Balke and Hamielec 
(6) is observed from Fig. 17b. 


The measured temperature history (Fig.l7a) for the NISO runs has been curve-fitted 


by: 

T{°C) = 37.S6 - 0.6S t + 0.27S6 for 0 < t ^ 8.3 min 

T{°C) = SO.O for t > 8.3 min 



im. No. 




T., KANPUR 

WJM 


( 12 ) 


where t is the time in minutes. These equations are used in the model of Ray et al. (23) to 
obtain theoretical predictions. Two sets of parameter values have been provided by Ray et 
al. (23). The individually optimized parameters (lOPs) have been obtained by curve fitting 
each of the isothermal experimental runs of Balke and Hamielec (6), and are available for 
every set of values of [I]o and T studied by them. These lead to excellent fits between model 
results and experimental data (6). On the other hand, best fit correlations (BFCs) have 
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(C) 

Fig. 17 Experimental results on the near isothermal (NISO) run 

(a) temperature history for the NISO a and b runs 

(b) conversion history for the NISO a and b runs 

(c) molecular weight vs. x for the NISOb run 

Theoretical predictions using the measured temperature history 
(Fig. 17a) and the BFCs in the model of Ray et al. (23) shown by 
solid curves. Isothermal results (BH) (6) also shown for comparision 
Dotted curve in Fig. 17b shows model results using the individually 
optimized paremeters (lOPs). [I]o=25.8 mol/m^, f°=0. 




55 


also been developed which relate the model parameters to the temperature in a functional 
form. These are more general since they can be used for values of [I]o and T which differ 
from those studied by Balke and Hamielec. These give curve fits which are not as good as 
those obtained using the lOPs, but are reasonably good. The lOPs and BFCs developed 
using isothermal data (6) are used without change for predicting nonisothermal results in 
our study. 

Equation 12. is used with BFCs to predict (solid curves in Fig.17 b and c) the 
monomer conversion and the number and weight average molecular weights (M^ and M^) 
as a function of time. The dashed curve in Fig. 17b gives the model (23) predictions for 
the near-isothermal run using the lOPs corresponding to T=50°C and [I]o=25.8 mol/m^. 
Good agreement is observed between our experimental data on monomer conversion and 
either of these predictions from the model (it may be noted that strictly speaking, we 
should not use the lOPs for cases where temperature is varying with time, as in the NI50 
runs). The experimental results on M„ are found (Fig. 17c) to lie close to the predicted 
values of using the BFCs, as expected. 

In the case of the NI70 run, the rise time (time taken to reach 70®C from room 
temperature ) is much larger, and hcis an important effect on the conversion history. Our 
results are found to deviate from those of Balke and Hamielec (6). Fig.lSb shows this 
deviation. The temperature history for this run is shown in Fig. 18a and has been curve- 
fitted by the following equations: 

T{°C) = 24.76 - 0.749 t + 0.5412 for 0 < t<: 10 min 

Ti^C) = 70.0 for t > 10 min (13) 

Model (23) results for conversion using both the sets of parameter values (lOP and 
BFC) with this temperature history have been generated and are shown in Fig.l8b. Fig.l8c 
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Molecular weight 





shows the experimental results and model predictions (using BFCs) for the molecular 
weights. The experimental values of M„ are observed to close to the predicted values of 
Mu;. Unfortunately, M„ data in the gel effect region could not be taken (this run is planned 
to be repeated) The x vs. t plot agrees reasonably well with model predictions. 

The next phase of the work was to obtain experimental data with step changes in 
temperature [step increase (SI) and step decrease (SD)] at different times and compare 
them with predictions from the model. Polymerization of MMA with the initiator (AIBN) 
loading of [I]o=25.8 mol/m® was carried out at 70°C (similar to NI70) for some time. After 
the system attained a steady temperature of 70°C, a step decrease of the temperature (SDl) 
at t=25 min (point a, Fig.lQa) to 50°C was effected using set point changes in the digital 
controller. The experimental results on conversion and along with the temperature 
history are shown in Figs.l9a-c. Two runs (SDla and SDlb) were made for this case to 
ensure reproducibility of the results with temperature transients. The temperature curve- 
fit equations are given by 

T{°C) = 24.7556 - 2.69 t + 1.28 f - 0.0629 f for 0 < t < 10.2 min 
T{°C) = 70.0 for 10.2 < t < 25 min 
T{°C) = 246 - 7 t for 25 < t < 28 min 

T{°C) = 50.0 fort >28 min (14) 

The model results are also shown in Figs. 19 b and c using the BFCs (23). It can 
be observed that the reaction slows down (due to lower rates at lower temperatures) sig- 
nificantly, and the gel effect is postponed compared to the NI50 run. A reasonably good 
agreement is observed between our experimental data and the model predictions. 
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In order to study the effect of a step decrease in the temperature at a later time, 
another polymerization run (SD2) was conducted under the same conditions as earlier but 
with the step decrease effected at t=45 min (point 6, Fig.20a). The curve-fit equations for 
the temperature history are given by 

Ti°C) = 18.069 - 0.9787 t + 0.934 - 0.0438 for 0 < t 9.8 min 

T{°C) = 70.0 for 9.8 < t < 45 min 

T{°C) — 246.713 — 4.297i for 45 < t 46.01 min 

T{°C) = 50.0 for t > 46.01 min (15) 

It can be seen that a delayed step decrease leads to a preponement of the gel effect as 
compared to the NI50 run (Fig. 20b). Again, a reasonably good agreement of our exper- 
imental results against model predictions using the BFCs (without any re-tuning of the 
parameters of the model) is observed in Figs.20 b and c. 

The effect of a step increase (SI) in temperature was also studied. Near isothermal 
polymerization was conducted at 50°C for some time and then, at a predecided time (t=120 
min, point c, Fig.21a) the temperature was increased to 70°C' by using a step change in 
the set point. Since the heating dynamics of the reactor system is a bit slow, the process 
took about 7 min to reach the desired set point. The experimental temperature history 
for this case was curve-fitted by 

T{°C) = 22.188 - 2.62392 t -t- 1.2611 - 0.046637 for 0 < t 7.13 min 

T(°C) = 50.0 for 7.13 < t < 120 min 

T{°C) = 6786.5 - 111.538 t + 0.4616 for 120 < t < 124.8 min 

T{°C) = 70.0 for t > 124.8 min (16) 

Model predictions using the best fit correlations and experimental x and data are 
shown in Figs.21 b and c respectively. Once again, reasonable accord between these is 
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Fig. 19 E.xperimental results on the step decrease (SDl a and b) runs. 

Theoretical predictions using the measured temperature history 
(Fig. 19a) and the BFCs in the model of Ray et al. (23) shown b\ 
solid curves. [l]o=25.8 mol/m^, f°=Q. 
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observed. 

Figs. 19b and 20b demonstrate experimentally the interesting result first predicted the- 
oretically by Ray et al. (23) that the gel effect can either get preponed or postponed when 
the temperature undergoes a step decrease, as compared to the case when the temperature 
is maintained constant (at the final value). The exact nature of the x vs. t plot depends 
on the viscosity (strictly speaking, the microviscosity) of the reaction mass at the time the 
step change is effected. Thus, if the polymer concentration and molecular weight at point 
a in Fig. 19b are such that the viscosity of the reaction mixture is lower than at point a’ 
on the NI50 curve, the gel effect will be postponed. Similarly, if the viscosity at point b 
(Fig.20b) is larger than at point b ’ on the NI50 curve, the gel ffect will be preponed. 

Our experimental results shown in Figs. 17-21 show an excellent accord with pre- 
dictions from the theoretical model of Ray et al. (23). It must be emphasized that the 
parameters (BFCs) in the theory were obtained using isothermal experimental data (6) 
taken on ampoule reactors, and were not re-tuned while predicting the present results us- 
ing experimental nonisothermal histories. This validates the model of Ray et al. (23) at 
least for nonisothermal reactor operations. Work is continuing to validate the theory for 
semibatch reactor operation, in which initiator /monomer are added to the reaction mass 
during polymerization. 
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CHAPTER 6 


CONCLUSIONS AND SUGGESTIONS FOR 

FUTURE WORK 


An experimental set-up to study the effect of temperature changes on bulk polymer- 
ization of MMA in a batch reactor has been assembled. A series of experiments under 
different temperature histories, at a fixed initiator loading, have been conducted. The ex- 
perimental results (conversion and average molecular weights) obtained show an excellent 
accord with predictions from the theoretical model of Ray et al. (23). This validates the 
model of Ray et al. (23) for nonisothermal reactor operations. It must be emphasized that 
the parameters in the model of Ray et al. were obtained using isothermal experimental 
data of Balke and Hamielec (6), and were not re-tuned for use in the present nonisother- 
mal study. These experimental results can be used to obtain better estimates of the model 
parameters. The present set-up can be used, with modifications, to operate the reactor 
under semi-batch conditions, in which initiator/monomer are added to the reaction mass 
during polymerization, and the results tested against the model of Ray et al. to validate 
the theory for semibatch reactor operations. 
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APPENDIX 1 


Table - A-1-1 

Conversion results for NI50a run 


Mass of monomer at t = 0 in reactor = 364.26 g 
Mass of initiator at t = 0 in reactor = 1.067 g 
Solvent fraction, f° = 0.0 ;[I]o = 25.8 mol/m^ 


No. 

Sampling 

mass of 

lOi+mass of 

t£; 2 +mass of 

mass of 

mass of 

monomer 


time 

empty 

benzene 

polymer 

polymer 

polymer 




sample 

bottle 

+ inhibitor 

sample 

sample 

precipitated 

conversion* 


(min) 

(g) 

^2 (g) 

^3(g) 

^‘'4(g) 

«^5(g) 


1 

25.00 

42.1818 

50.8478 

51.0720 

0.2242 

0.0029 

0.0131 

2 

45.00 

63.1956 

71.8840 

73.3340 

1.4500 

0.0989 

0.0682 

3 

60.00 

63.8669 

72.5351 

76.6868 

4.1517 

0.2960 

0.0713 

4 

80.00 

54.4975 

63.1141 

71.5454 

8.4313 

0.8364 

0.0992 

5 

100.00 

51.8205 

60.4698 

63.6862 

3.2164 

0.4310 

0.1340 

6 

125.00 

49.8981 

55.5580 

61.6065 

6.0485 

0.9133 

0.1510 

7 

140.00 

55.0350 

63.6920 

65.1961 

1.5041 

0.2602 

0.1730 

8 

150.00 

50.5054 

59.1649 

59.5500 

0.3851 

0.0728 

0.1890 

9 

180.00 

42.2414 

50.9066 

51.1360 

0.2294 

0.0548 

0.2390 

10 

186.00 

61.7470 

70.4100 

73.6900 

3.2800 

0.8561 

0.2610 

11 

190.00 

62.2074 

70.8570 

75.9600 

5.1030 

1.6483 

0.3230 

12 

193.00 

61.6740 

70.3488 

74.6400 

4.2912 

1.4333 

0.3340 

13 

200.00 

58.7476 

67.4358 

76.6340 

9.1982 

3.1090 

0.3380 

14 

220.00 

64.6650 

73.3320 

78.5892 

5.2572 

2.3710 

0.4510 

15 

225.00 

64.9650 

73.6532 

79.6530 

5.9998 

2.8919 

0.4820 

16 

230.00 

61.9486 

70.6368 

78.5380 

7.9012 

4.0533 

0.5130 

17 

238.00 

62.4820 

71.1744 

78.2360 

7.0616 

5.0279 

0.7120 

18 

260.00 

61.5632 

70.2514 

75.2130 

4.9616 

4.2322 

0.8530 


conversion x 

Wi, 
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Table - A-1-2 

Conversion results for NI50b run 


Mass of monomer at t = 0 in reactor = 541.72 g 
Mass of initiator at t = 0 in reactor = 2.394 g 
Solvent fraction, f° = 0.0 ;[I]o = 25.8 mol/m^ 



Sampling 

mass of 

tui+mass of 

u; 2 +mass of 

mass of 

mass of 

monomer 

■ 

time 

empty 

benzene 

polymer 

polymer 

polymer 


■ 


sample 

bottle 

+ inhibitor 

sample 

sample 

precipitated 

conversion* 

■ 

(min) 

(g) 

W2 (g) 

^«3(g) 

Mg) 

«^5(g) 


1 

25.00 

42.2474 

50.8801 

58.9550 

8.0749 

0.1130 

0.0140 

2 

50.00 

56.9174 

65.5693 

67.5825 

2.0132 

0.1278 

0.0635 

3 

75.00 

51.8283 

60.4727 

65.0394 

4.5667 

0.3745 

0.0820 

4 

100.00 

48.6316 

57.2556 

60.9876 

3.7320 

0.5038 

0.1350 

5 

150.00 

61.6728 

70.3271 

71.8861 

1.5590 

0.2853 

0.1830 

6 

175.00 

39.6350 

48.2803 

52.3484 . 

4.0681 

0.8868 

0.2180 

7 

200.00 

48.9352 

57.5823 

59.1000 

1.5177 

0.4978 

0.3280 

8 

210.00 

58.8449 

67.4713 

67.7429 

0.3209 

0.1107 

0.3450 

9 

212.00 

64.1432 

72.7765 

73.8853 

1.1088 

0.3845 

0.3468 

10 

220.00 

51.0863 

59.0520 

60.6854 

1.6334 

0.6534 

0.4000 

11 

223.00 

49.2256 

57.8671 

58.2238 

0.3567 

0.1569 

0.4400 

12 

230.00 

62.2216 

70.8113 

70.9363 

0.1250 

0.0585 

0.4680 

13 

235.00 

49.3197 

57.9417 

58.9563 

1.0146 

0.4870 ; 

0.4800 

14 

240.00 

51.1434 

59.7930 

60.1611 

0.3681 

0.1914 

0.5200 

15 

242.00 

65.4352 

74.7545 

81.6118 

6.8573 

5.8973 

0.8600 


conversion x 

Wi 
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Table - A-1-3 

Conversion results for NI70 run 


Mass of monomer at i = 0 in reactor = 373.6 g 
Mass of initiator at < = 0 in reactor = 1 .6924 g 
Solvent fraction, f° = 0.0 ;[I]o = 25.8 mol/mr 



Sampling 

mass of 

lOi+mass of 

ti; 2 +mass of 

mass of 

mass of 

monomer 

■ 

time 

empty 

benzene 

polymer 

polymer 

polymer 


■ 


sample 

bottle 

+ inhibitor 

sample 

sample 

precipitated 

conversion* 

■ 

(min) 

wi (g) 

W2 (g) 

«^3(g) 

^"4(g) 

^"5(g) 


1 

10.00 

41.4161 

50.1764 

"68.2180 

18.0416 

0.1967 

0.0109 

2 

20.00 

41.5387 

50.2505 

55.3020 

5.0515 

0.4551 

0.0901 

3 

25.00 

56.7561 

65.4333 

73.5700 

8.1367 

0.9064 

0.1114 

4 

30.00 

49.2877 

57.9870 

64.5700 

6.5830 

0.9749 

0.1481 

5 

35.00 

49.2972 

56.9795 

64.1130 

7.1335 

1.2491 

0.1751 

6 

40.00 

49.2976 

57.9609 

64.5460 

6.5851 

1.4171 

0.2152 

7 

45.00 

62.5471 

71.2286 

75.5664 

4.3378 

1.0090 

0.2326 

8 

50.00 

58.8406 

67.5330 

69.4391 

1.9061 

0.5305 

0.2783 

9 

53.00 

62.5241 

71.2989 

74.1855 

2.8866 

1.1087 

0.3841 

10 

54.00 

50.5061 

59.1823 

61.8772 

2.6949 

1.1752 

0.4361 

11 

56.00 

49.2356 

57.9206 

59.8110 

1.8904 

0.9760 

0.5163 

12 

58.00 

59.3727 

68.0559 

69.0857 

1.0298 

0.5851 

0.5682 

13 

60.00 

49.8990 

58.5703 

59.8832 

1.3129 

0.7894 

0.6013 

14 

80.00 

56.3214 

65.0122 

66.8011 

1.7889 

1.5664 

0.8756 


conversion a: 

tt'4 
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Table - A-1-4 

Conversion results for SDla run 


Mass of monomer at t = 0 in reactor = 373.6 g 
Mass of initiator at t = 0 in reactor = 1.6512 g 


Solvent fraction, f° = 0.0 ;[I]o = 25.8 mol/m' 


No. 

Sampling 

mass of 

Wi+mass of 

i/; 2 +mass of 

mass of 

mass of 



time 

empty 

benzene 

polymer 

polymer 

polymer 




sample 

bottle 

+ inhibitor 

sample 

sample 

precipitated 



(min) 

(g) 

^^^2 (g) 

U^3(g) 

^^4(g) 

“’slg) 


1 

5.00 

42.1818 

50.8881 

52.7634 

1.8753 

0.0034 

0.0018 

2 

10.00 

63.8669 

72.5730 

74.5840 

2.0110 

0.0402 

0.0200 

3 

20.00 

63.1956 

71.9019 

74.7360 

2.8341 

0.2281 

0.0805 

4 

30.00 

54.4975 

63.2038 

65.3850 

2.1812 

0.3300 

0.1513 

5 

40.00 

51.8205 

60.5268 

61.6380 

1.1112 

0.1948 

0.1753 

6 

50.00 

49.8981 

58.6044 

61.3230 

2.7186 

0.4385 

0.1613 

7 

60.00 

55.0350 

63.7413 

64.9876 

1.2463 

0.2297 

0.1843 

8 

90.00 

42.2414 

50.9477 

54.8360 

3.8883 

0.7971 

0.2050 

9 

120.00 

50.5050 

59.2113 

63.8320 

4.6207 

1.1062 

0.2394 

10 

150.00 

49.3483 

58.0546 

59.9830 

1.9284 

0.5288 

0.2742 

11 

170.00 

62.2074 

70.9137 

74.9137 

4.0000 

1.2212 

0.3053 

12 

180.00 

61.6740 

70.3803 

72.1380 

1.7577 

0.5285 

0.3007 

13 

190.00 

63.5631 

72.2694 

74.2830 

2.0136 

0.6393 

0.3175 

14 

210.00 

58.7863 

67.4926 

69.9530 

2.4604 

0.8611 

0.3500 

15 

218.00 

59.8363 

68.5966 

70.2340 

1.6374 

0.5777 

0.3528 

16 

220.00 

63.5463 

72.3066 

73.9870 

1.6804 

0.5975 

0.3556 

17 

225.00 

63.4320 

72.1383 

75.3840 

3.2457 

1.2282 

0.3784 

18 

230.00 

59.7634 

68.4697 

70.4830 

2.0133 

.7959 

.3953 


conversion x 

tUi 
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Table - A-1-5 

Conversion results for SDlb run 



5.00 

10.00 
20.00 

25.00 

40.00 

90.00 

150.00 

170.00 

210.00 
225.00 

250.00 

270.00 

275.00 

280.00 


Mass of monomer at 
Mass of initiator at t 
Solvent fraction, /° = 

mass of mi+mass of 
empty benzene 
sample + inhibitor 


t = 0 in reactor 
= 0 in reactor = 
= 0.0 ;[IL = 25, 


empty 

sample 

bottle 

(g) 

42.1207 

63.1415 

63.8127 

49.8441 

51.7781 

42.1874 

49.2853 

62.2238 

54.6779 

61.2253 

58.9251 

61.1047 
58.9251 

60.1047 


(g) 

50.8810 

71.9018 

72.5730 

58.6044 

60.5384 

50.9477 

58.0456 

70.9841 

63.4382 

69.9856 

67.6854 

69.8650 
67.6854 

68.8650 



= 513.7 g 
: 2.2704 g 

8 mol/m_^ 

mass of 
polymer 
precipitated 


52.3840 

74.5320 
75.8360 
63.8430 
62.2386 
54.9860 
63.2840 
73.8620 

68.5320 
71.9560 
70.8320 
72.9350 
70.8342 
71.9216 
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Table - A-1-6 

Conversion results for SD2 run 


Mclss of monomer at i = 0 in reactor = 373.6 g 
Mass of initiator at t = 0 in reactor = 1.6512 g 
Solvent fraction, f° = 0.0 ;[I]o = 25.8 mol/m^ 


No. 

Sampling 

mass of 

lUi+mass of 

ti; 2 +nnass of 


mass of 

monomer 


time 

empty 

benzene 

polymer 

polymer 

polymer 




sample 

bottle 

+ inhibitor 

sample 


precipitated 

conversion* 


(min) 

(g) 

'^2 (g) 


«^4(g) 

^^5(g) 


1 

10.00 

63.2631 

72.0234 

75.5857 

3.5623 

0.0363 

0.0102 

2 

20.00 

56.9517 

65.7120 

67.6157 

1.9037 

0.2005 

0.1053 

3 

30.00 

41.5674 

50.3277 

52.2001 

1.8724 

0.2831 

0.1512 

4 

40.00 

51.2778 

60.0381 

65.5886 

5.5505 

1.1634 

0.2096 

5 

50.00 

49.4287 

58.1890 

61.0580 

2.8690 

0.9313 

0.3246 

6 

60.00 

51.1781 

59.9384 

61.0450 

1.1066 

0.3714 

0.3356 

7 

75.00 

65.2048 

73.9651 

76.0670 

2.1019 

0.7085 

0.3371 

8 

120.00 

62.4303 

71.1906 

71.5672 

0 .3766 

0.1368 

0.3632 

9 

150.00 

62.1590 

70.9193 

71.3410 

0.4217 

0.1692 

0.4013 

10 

210.00 

49.1239 

57.8842 

59.0101 

1.1259 

0.5813 

0.5163 

11 

215.00 

52.2817 

61.0420 

64.9915 

3.9495 

2.1375 

0.5412 

12 

218.00 

40.7731 

49.5334 

53.0078 

3.4744 

1.9568 

0.5632 

13 

222.00 

15.5167 

24.2770 

28.2581 

3.9811 

3.2366 

0.8130 


conversion x 

Wi 
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Table - A-1-7 

Conversion results for S% run 


Mass of monomer at t = 0 in reactor = 467.0 g 
Mass of initiator at t = 0 in reactor = 2.064 g 
Solvent fraction, = 0.0 ;[I]o = 25.8 mol/m^ 


No. 

Sampling 

mass of 

wi+mass of 

ti) 2 +mass of 

mass of 

mass of 

monomer 


time 

empty 

benzene 

polymer 


polymer 




sample 

bottle 

+ inhibitor 

sample 

sample 

precipitated 

conversion"* 


(min) 

(g) 

'^2 (g) 


^^4(g) 

wsig) 


1 



59.7975 

63.9566 

4.1591 

0.0524 

0.0126 

2 



69.6237 


7.7167 

0.1173 

0.0152 

3 


52.4835 

61.2438 

65.1651 

3.9213 

0.2012 

0.0513 

4 




78.6636 

8.6363 

0.6598 

0.0764 

5 



59.7465 

73.7833 

14;0368 

1.1819 

0.0842 

6 


62.4883 

71.2486 

84.8275 

13.5789 

1.7381 

0.1280 

7 



65.0397 


9.0583 

1.2464 

0.1376 

8 

1 

115.00 

49.2732 

58.0335 

66.6431 

8.6096 

1.3362 

0.1552 

9 

125.00 

63.4426 

72.2029 

83.7706 

11.5677 

3.1383 

0.2713 



42.2717 

51.0320 


6.0739 

2.6276 

0.4326 

11 

132.00 

41.3518 

50.1121 

55.3437 

5.2316 

2.4243 

0.4634 

12 

135.00 

43.3216 

52.0512 


8.0853 

4.4193 

0.5132 


conversion x 

Wi 
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Molecular weight results for Nl50b run 


No. 

Polymer 

concentration** 

C. (g/dl) 

Cl 

Efflux t 

C 2 

imes (s) 

C3 

C4 


X 

time 

min 

1 

0.75 

341.93 

289 

235.47 

215.7 

8.162e+05 

0.135 

100 

2 

1.0 

390.5 

341.63 

270 

234.6 

8.1064e+05 

0.183 

150 

3 

1.25 

373.74 

380.92 

240.0 

203.65 

2.963996-f-05 

0.218 

175 

4 

0.745 

341.0 

289.0 

235.77 

215.7 

7.70112e+05 

0.328 

200 

5 

0.745 

341.93 

289.06 

235.473 

215.7 

8.319e+05 

0.345 

210 

6 

0.832 

410.34 

340.76 

260.45 

200.5 

4.927e+05 

0.3468 

212 

7 

0.75 

548.458 

405.001 

289.684 

242.581 

9.38955e+05 

0.4 

220 

8 

0.84 

698.261 

496.742 

337.126 

273.033 

1. 06764 9e+06 

0.44 

223 

9 

0.8151 

405.65 

326.85 

282.46 

— 

12.9e+05 

0.468 

230 

10 

1.000 

231.61 

207.98 

184.72 

175.97 

2.40e+05 

0.48 

235 

11 

1.0 

400.34 

310.23 

230.65 

200.76 

3.14022e+05 

0.52 

240 

12 

0.819 

491.73 

367.93 

269.12 

228.23 

7.30049e+05 

0.86 

242 


efflux time for pure solvent (benzene), to = 143.6s 


Molecular weight results for NI70 run 


No. 

Polymer 


Efflux times (s) 


M„ 

X 

time 


concentration** 
Cl (g/dl) 


C2 

C 3 

C4 



min 

1 

0.3045 

189.185 

177.655 

164.145 

159.11 

4.62385e+5 

0.0183 

12.56 

2 

0.878 

256.66 

222.73 

193.565 

179.05 

3.13932e+5 

0.11498 

25 

3 

1.3215 

297.64 

251.165 

210.5 

191.385 

2.54268e+5 

0.1315 

30 

4 

0.9923 

244.39 

216.68 

191.225 

173.595 

2.32986e+5 

0.2049 

35 

5 

0.922 

230.04 

211.49 

281.75 

171.255 

2.21975e+5 

0.2475 

40 

6 

0.8376 

236.305 

209.485 


173.475 

2.85087e+5 

0.2481 

42 

7 

0.7563 

234.335 

209.715 

185.785 

176.395 

3.65839e+5 

0.2522 

46.5 

8 

0.8033 

228.265 

203.745 

182.265 

171.385 

2.75924e+5 

0.2648 

48.5 

9 ^ 

0.208 

158.9 

154.25 

149.455 

147.525 

2.57805e+5 

0.2940 

50 


0.8803 

230.1 

204.97 

180.775 

169.475 

2.15449e+5 

0.3195 

53 


efflux time for pure solvent (benzene), to = 147.08s 

C2 = 15 Cl/20 
C 3 = IOC2/I5 
C4 = I5C3/2O 
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Molecular weight results for SDla run 


No. 

Polymer 
concentration** 
cr (g/dl) 

Cl 

Efflux ti 

C2 

mes (s) 

C3 

C4 


X 

time 

min 

1 

0.745 

341.43 

289.56 

235.44 

215.77 

8.85762e+05 

0.020 

10 

2 

1.0 

232.205 

207.53 

185.21 

176.68 

2.3967e+05 

0.0805 

20 

3 

1.0 

400.64 

312.23 

230.42 

201.02 

3.14132e+05 

0.1513 

30 

4 

1.0 

400.34 

310.23 

230.65 

200.76 

3.14022e+05 

0.1754 

40 

5 

1.0 

410.43 

318.54 

240.83 

200.32 

3.05443e+05 

0.2050 

90 

6 

1.0 

413.4 

320.5 

243.4 

210.3 

4.30934e+05 

0.3175 

190 

7 

1.0 

421.56 

325.52 

246.53 

214.34 

4.66681e+05 

0.3500 

210 

8 

0.832 

410.23 

341.23 

260.34 

200.78 

4.95676e+05 

0.3528 

218 

9 

1.0 

468.54 

354.43 

262.53 

224.87 

5.44117e+05 

0.3556 

220 

10 

0.818 

490.83 

368.23 

269.12 

228.23 

7.27149e+05 

0.3784 

225 

11 

0.745 

341.8 

288.0 

235.4 

215.4 

8.31302e+05 

0.3953 

230 


efflux time for pure solvent (benzene), to = 143.58s 


Molecular weight results for SDlb run 


No. 

Polymer 


Efflux times (s) 



X 

time 


concentration** 
Cl (g/dl) 

Cl 

C2 

C3 

C4 



min 

1 

1.00 

232.59 

208.00 

184.62 

176.97 

2.3998e-l-05 

0.0793 

20 

2 

1.00 

232.61 

207.98 

185.62 

176.97 

2.3987e+05 

0.1037 

25 

3 

1.0 

400.64 

312.22 

230.39 

201.00 

3.14e+05 

0.1613 

40 

4 

1.0 

410.39 

318.53 

240.79 

200.30 

3.056e+05 

0.1973 

90 

5 

1.0 

410.43 

318.54 

240.83 

200.32 

3.05443e+05 

0.2831 

150 

6 

1.0 

413.4 

320.5 

243.38 

210.29 

4.308e+05 

0.3011 

170 

7 

0.830 

410.23 

341.23 

260.30 

200.72 

4.95e+05 

0.3526 

210 

8 

1.0 

468.54 

354.43 

262.53 

224.77 

5.45e+05 

0.4367 

270 

9 

0.819 

491.53 

367.91 

270.02 

229.13 

7.31e+05 

0.5213 

275 

10 

0.744 

342.0 

289.0 

235.4 

215.7 

8.32e+05 

0.5626 

280 


efflux time for pure solvent (benzene), to = 141.63s 
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Molecular weight results for SD2 run 


No. 

Polymer 

concentration** 

c. (g/dl) 

Cl 

Efflux 1 

C2 

,imes (s) 

C3 

C4 


X 

time 

min 

1 

1.0043 


207.9 

184.30 

173.56 

2.0745e+05 

0.1512 

30 

2 

1.005 

-232.61 

207.98 

185.62 

176.97 

2.3887e+05 

0.2096 

40 

3 

1.0 


318.54 

240.83 

200.32 

3.05443e+05 

0.3246 

50 

4 

1.0 

421.56 

325.52 

246.53 

214.34 

4.66681e+05 

0.371 

75 

5 

0.832 


341.23 

260.34 

200.78 

4.95e+05 

0.3632 

' 120 

6 

1.0 


310.23 

230.65 

200.76 

3.139e+05 

0.4013 

150 

7 

0.819 

491.73 

367.93 

269.12 

228.23 

7.29e+05 

0.5412 

215 

8 

0.745 

341.93 

289.06 

235.473 

215.7 

8.319e+05 

0.5632 

218 

9 

0.745 

341.43 

289.56 

235.44 

215.77 

8.85e+05 

0.813 

222 


efflux time for pure solvent (benzene), to = 140.63s 


Molecular weight results for SI run 


No. 

Polymer 
concentration** 
Cl (g/dl) 

Cl 

Efflux t 

C2 

imes (s) 

C3 

C4 


X 

time 

min 

1 

1.0 

390.49 

341.629 

269.99 

234.6 

8.1056e+05 

0.0152 

30 

2 

0.745 

341.43 

289.56 

235.44 

215.77 

8.857e+05 

0.0513 

45 

3 

0.745 

341.93 

289.06 

235.473 

215.7 

8.318e+05 

0.0764 

60 

4 

0.745 

341.0 

289.0 

235.77 

215.7 

7.6996e+05 

0.0842 

75 

5 

0.75 

341.93 

289 

235.47 

215.7 

8.162e+05 

0.128 

90 

6 

0.745 

341.0 

288.99 

235.6 

214.7 

7.727e+05 

0.1376 

105 

7 

0.745 

341.0 

289.0 

235.77 

215.7 

7.701 12e+05 

0.1552 

115 

8 

0.75 

548.45 

405.0 

289.64 

242.58 

9.3955e+05 

0.2713 

125 

9 

0.744 

342.0 

289.0 

235.4 

215.7 

8.32e+05 

0.4634 

132 

10 

0.84 

697.251 

496.740 

337.123 

273.003 

1.075e+06 

0.5132 

135 


efflux time for pure solvent (benzene), to = 143.62s 
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APPENDIX 2 


FLOW CHflHT FOR THE CONTROL PROGRAM 


CONTROL I 


START I 




N.B I Proctdur* ar» srlvtn In quotations 








8 
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C 


HEAD_I H_COHTROL 
DISPLAy_COMTROL 












'REflD_lH_BTEP' I 
'D16PLfty_6TEP' 1 





REPEAT 




























APPENDIX 3 


DETERMINATION OF M„ 


Molecular weight Mt, is obtained from the efflux times as described below (36-38): 

Let Cl be the concentration of the original solution prepared (all concentrations in 
g polymer/dl solution). Dilution of 15 ml of this solution with 5 ml of benzene gives 
C 2 = 15(71/20. 10 ml of this solution is diluted with 5 ml of benzene, to give Cz — 10(72/15. 
Dilution of this 15 ml solution with 5 ml of benzene results in C 4 = 15(73/20. Let the 
corresponding efflux times be t,-, i=l,....,4 (averages of 4 readings which are within about 
0.5s). Let to be the efflux time for the pure solvent. 

The relative viscosity, specific viscosity, reduced viscosity and inherent viscosity for 
the solution are defined as follows: 


Vr,i 

= U/to 

( 1 ) 

Vsp,i 

^ Vri 1 

( 2 ) 

Vred^i 


( 3 ) 

Vinhyi 

= ln{rir,i/Ci) 

( 4 ) 


The intrinsic viscosity (which is related to the viscosity average molecular weight) is de- 
fined as the reduced viscosity in the limit of vanishing concentration. Mathematically, 

[»/] = Jjm Tired (5) 

At low concentrations (below about 0.5 g/dl), plots of reduced viscosity vs concentration 
and of inherent viscosity vs concentration are linear, following the established equations: 

Vred = [rj] + k[jifc 
rjinh = [»?] - k[rif c 


( 6 ) 

( 7 ) 
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Therefore, it can be observed that botli these linear functions have a common intercept 
at the intrinsic viscosity, [t;]. A double extrapolation routine which fits Equations (A3. 6, 
A3. 7) in a least squares sense with a common intercept was used to obtain [r}]. A listing of 
the program is included at the end of this appendix. A specimen plot of the determination 
of [//] .is shown in Fig A3.L 

The M„ is calculated from the Mark-Houwink equation, 

W = KMt (S) 

The values of the constants, K and a, were taken from the Polymer Handbook (39) for 
the PMMA -Benzene system at 30°C. The last run samples were analyzed at 32'’C (in a 
water bath), as the room temperature had gone above 30°C. Still the same constants were 
employed. This is not expected to lead to much little error as these constants, K and a, 
are not sensitive functions of temperature (39). The following values of K and a are used 
for the calculations; 

K = 5.2 X 10"® dl/g 


a 


0.76 




Fig. A3.1 Double extrapolation plot to estimate the intrinsic viscosity 
(of sample 4, NI50b run). = 1.618 dl/g, Mv= 8.162 x 10® 


